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Abstract 

An elegant and more precise (Denominator) formula for the 3-loop perturbative 
QCD coupling is discussed. It improves the common expression (e.g., canonized by 
PDG) in few GeV region. On its base, we propose simple analytic Model for ghost- 
free QCD running couplings and their effective powers within the Analytic Per- 
turbation Theory, in both the space-like (Euclidean) and time-like (Minkowskian) 
regions, very accurate in the range above 1 GeV. 

Effectiveness of the new Model is illustrated by the example of T(1S) decay 
where the standard analysis gives as{M~£) = 0.170ib0.004 value that is inconsistent 
with the bulk of data for q^. Instead, we obtain af°'^{Mr) = 0.185 ± 0.005 that 
corresponds to a^^°'^{Mz) = 0.120 ± 0.002 that is close to the world average. 

1 Introduction 

This text contains essence of three topics related to the QCD coupling at low energy: 

1. Presenting of QCD perturbative coupling as(/i) by Denominator representation in- 
stead of "canonical" PDG-like bulky expression. 

2. Particular "Analytic-Perturbation-Theory" (APT) ghost-free model for as{fi) in the 
low-energy region. 

3. APT Model analysis of Upsilonium decay. 

It is based mainly upon recent paper which contain more detailed exposition. 

1.1 Notation 

We use the Bethke[2j notation (different from the PDG one) 

^=/?(«,) = -/?o«^-A«'-/52«'+ ... , L = \n{x/A'), (1) 

for beta-function coefficients [3 (a) = —(3q (1 + 6i a + 62 + . . . ) bk = Pk/ Po , 

R( A 33 - 2/ , , , 153 - 19/ ,MS, , 2857 - (5033/9)n^ + (325/27)n^^ 
^°^"^) = ^^' ^^^"^) = 2vr(33-2/)= ^'^^ = 32.^(11 - (2/3)n,) 
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Here, numerically, all (3k , and B = /Sf/jSo are of an order of unity 

/3o(4 T 1) = 0.6631 ± 0.0530 ; /3i(4 =F 1) = 0.3251 ± 0.0802 ; /32(4 =F 1) = 0.205l2\^i°2 ; 
6i(4 ± 1) = 0.4902;0:08f7 . ^^(4 ^ ^ 0.309;0:}i9 ; B{A T 1) = 0.7392t[|:0^09 . 

2 Perturbative QCD coupling 

2.1 Denominator representation for 3-loop as{n) 
Instead of canonical ( Bethke[2], PDG[3]) 1/L expanded expression 

= ^ - 1^ + L - In L - 1) + 62] - 

[6?(ln3L-fln2L-21nL + i) +36i62lnL-|] +... {9.5pdg) 
we argue to use the "Denominator representation" for a;s(yu) 

^^"'"'^(L) = , \ -—; L = ln^ (2) 



Canonical form can be obtained from it by use of casual expansion parameter €2 

{3,iter) T \ _ 1 I J I ^2 _ c„ _L £4 



^siL) = ^ ( 1 + —r + 4-e3 + 



1 + (62 + 6262 + el) - (£3 + ^^ 



2, 



-1 



B InL 
= — « 1 



as well as the higher loops ones: e^+i = jj^^ > 1) • For instance, at n/ = 3, 4 

.^ r^-,^ ^-..ilni .^ [0.79^-0.74] ,^ /0.90-f0.84\^ /1. 2 

€2(3-4) = [0.79-0.74] — ; £2(3^4) = ^ i; £3(3-4) = f j; £4 = f — 

Then, ea^ ~ £3 at InL* ~ 1.13 : L* ~ 3.3, that is at ~ 5.0 A, Q ~ 1.6 ^2GeV. 

There, 62 ^ 0.25 , £3 ~ 0.07 . Relative {Q.^pdg) error ~ €2^ ~ 0.02 , while error of (2) 
~ £4/2 ~ 0.03 Ghost singularity at L = in ([2]) is much weaker than, in {^.^pdg)- 
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2.2 The beta function arbitrariness 



In turn, Denominator representation ([2]) is an iterative approximate solution of RG dif- 
ferential eq.(l) quadrature which, in the two-loop (NLO) case, is 

(^oL = - 2(,,u \ = — + ^1 In — T 3) 

J a^(l + Oia) as 1 + Oi 

- a transcendental relation that can be resolved [1] in terms of special Lambert function. 
Usually, instead, one solves this relation iteratively with the result 

' ^ ' PoL + h HB + l/a^'^Po) f^oL + B(3o ln(L + B) ^ ^ 

with its famous log-of-log dependence. 

Meanwhile, if one transforms 2-loop beta function a la Pade (32 — > (32{ct) ~ — 
the result of integration and iteration will be different 



l-bi a ' 



(7(1 1 1 

(l_6^a) = - + 6iln(/3oa.), c^T'^'^'iL) = —-——- . (5) 
as Po + bi m L 

Thus, the inner self- consistency condition of the 2-loop (NLO) approximation looks like 

Aa(2) a(2P)_c,{2) 



a(2) 



^2 « 1 ■ (6) 



Due to this, the intrinsic accuracy estimate for the two-loop QCD coupling a*^^) at 
1.5 - 2 GeV region is about 7 %, that is (in the MS scheme) equal to the value of 
three-loop contribution ~ £3 . At the same time, at the b-quark mass ~ 4.5 GeV one has 
el - 0.015, 62 ~ 0.21, ~ 0.01 

This means that in the nj = 4,5 domain one can use only the first line of the 3-loop 
eq.(9.5pDG) and second expression (5) instead of (2) with the same accuracy. 

Numerically, at 4-flavor region for A4 ~ 300 MeV and q = 1.5; 3 GeV one has 62{q = 
1.5 GeV) ~ 0.053, 63^(1. 5) ~ 0.020,62(3) ~ 0.025, e2^(3) ~ 0.020 and at 5-fiavors 
for A5 ~ 225 MeV ; g ~ 4.5 GeV - ^2(4.5) ~ 0.012 , £^(4.5) ~ 0.030 . 

Pragmatically, this means that, within the 2 per cent limit of accuracy, one can equally 
use simple two-term denominator eq.(5) up to 3 GeV scale. However, in the 1/L-expanded 
form one needs to keep 3 terms (of expansion in powers of €2) in = 5 region and 5 
terms in the Uf = 4 one. 

As it can be shown, the inner consistency of the 3-loop iterative approximations is 
controlled by the same parameters 62 and €2 ■ At the same time, the measure for impor- 
tance of direct 3-loop (NNLO) term = b2/{Po should be compared with 62 ■ It turns 
out that they are close (within 40 per cent) to each other. 
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3 Analytic Perturbation Theory 

3.1 Outline of Analytic Perturbation Theory 

Remind first, that the cornerstones of APT are the analyticity of coupling functions and 
compatibility with linear integral transformations. For the fresh reviews of APT see [5l[6]. 

Here follows compendium of main definitions. The most elegant APT formulation is 
based on the set of spectral functions {pj(cr)} defined as 

p,iz) = lmi[asi-z)]'). (7) 

The first of them, pi = p(cr) is just the Kallen-Lehmann spectral density for the Euclidean 
APT coupling. Then, higher Euclidean ( "analyticized kth power of coupling in the Eu- 
clidean domain" ) and Minkowskian ( "effective kth power of coupling in the Minkowskian 
domain" ) APT functions will be 

A,{Q') = A[a';] = - n ■ Ms) = ffi[a.] = - H -Pk{a) . (8) 

Jo + Q J, CT 

They are related by integral transformation and satisfy differential relations 

which can be used for iterative definitions. 

3.2 Properties of the APT functions 

For the one-loop case, the APT formulae are simple and elegant. Starting with the per- 

turbative RG-improved QCD coupling ai'^\Q'^) = l/(/?oO > with the help of ([7]), (jH]) one 
arrives at the ghost-free effective Euclidean and MinkowskiarQ 

APT couplings. Higher functions Ai , Stj can be defined via recursive relations ([9]) with 
only one term in the r.h.s. For instance, Minkowskian functions are 

^2'^(^) = i72^' ^3(^) = -^ ,,,^ ,,, , 2t4(L) = f ~ 3 . (11) 



^Here, we change arguments of the APT functions: I ~ hi(Q^/A^) and s L ~ ln(s/A-^) . 
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Higher-loop case. The two-loop expressions are more complicated. Here, exact QCD 
coupling Us can be expressed explicitly in terms of a special Lambert function W defined 
as a solution of the transcendental equation W{z) = z . This expression yields rather 
complicated formulae for a^; = Ai and aq = 2li in terms of the main branch W-i . 

At the three-loop case, one meets further complications. Here, only for Fade approxi- 
mated beta-function, exact solution can be expressed [4j in terms of the Lambert function. 
Such expression are not comfortable enough for practical use. 

The devised scheme with due account for matching [7], known as "global APT", has 
been studied numerically by Magradze and Kourashev at the two- and three-loop level. 
They calculated numerical tables for the first three functions 2tfcl, 2, 3 and ^1,2,3 at three 
values of A"/=3 = 350, 400, 450 MeV in the interval 1 GeV < ^/s,Q < 100 GeV' [l [8], and 
2li,2 , Ai,2 in the interval 0.1 GeV <^/s,Q <3 GeV [9J. 

The APT functions obey important properties valid in the higher-loop case: 

• In the Euclidean and Minkowskian domains, QCD couplings a^lQ^) = Ai, (Q^) , 
«m(s) = 2li and their "effective powers" Ak{Q'^) 2lfc(s) , are different functions 
related by integral operations A[ ] and ]R[ ] explicitly defined in eqs. (E]). Higher 
functions like (8), (9), are not equal to powers of the first ones (ITOl) . 

• The APT functions Ak{Q^) 2t/c(s) differ of common expansion functions (tts)'^ in 
the low energy region, where they are regular with finite aE^O) = ^^(O) = l/Po 
or zero limits. Unphysical singularities are absent with no additional parameters 
introduced. This behavior provides high stability with respect to change of renor- 
malisation scheme [10]. In the UV limit, all APT functions tends to their usual 
counterparts (ag)^ . 

• Expansion of an observable in coupling powers (^^(Q^))" for the Euclidean or in 
{as{s))^ for the Minkowskian case is substituted by nonpower expansion in sets 
{Ak} , or {2tfc} respectively. The latter expansions exhibit a faster convergence. 

The APT re-examination of various processes has been performed in number of pa- 
pers pjj - [IA\. In particular, paper pjj contained the first attempt in revising the T 
decay. Simple estimates for the influence of tt^ — terms upon some observables were per- 
formed in [5]. The r decay was re-examined in [12]. Pion form factor was studied [T3| [Ti] 
within the APT techniques. One more extraction of APT coupling at Q ~ 100 — 400 MeV 
was made [15] by the Milano group from mass spectrum analysis of ground and first excited 
quarkonium states. 
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4 Simple Model for 3-Loop APT Functions 



4.1 "One-Loop-Like" Model 

Our aim is to construct simple and accurate enough (for practical use) analytic approx- 
imations for two sets of functions 21^ and Ak , k = 1,2,3 . To reduce number of fitting 
parameters, one should better provide the applicability of the recurrent relations. To 
this goal, we use one-loop APT expressions, eqs. lfTOll - ( fTTl) . with modified logarithmic 
arguments 

Ar'il) = Ai'\Q ; 2ir'(^) = , (12) 

and being some "two-loop RG times". Model functions (fT2l) are related by the 
"one-loop-type" recursive relations 

1 J A mod 1 J A mod j; i r/9i'"'"^ 

Amod _ ^ "-^n _ ^ "-^n ^mod _ ^ ""^n 



A simple expression for was taken from [TB], where a plain approximation for the 
two-loop effective log I2 = / + 6 In ^/W^^An^ , with b defined in Section 1 was used. This 
approximation combined reasonable accuracy in the low-energy range with the absence of 
singularities for aE ■ We extend this approach to higher functions in both the Euclidean 
and Minkowskian domains and change square root in "effective logs" -^2(0) and hia) '■ 
VWT4^ VP + an^ with a , an adjustable parameter. It comes out from thorough 
numerical analysis that optimal value of the new parameter is a ~ 2 , while effective 
boundaries between the flavor regions have to be chosen on quark masses nic = 1.3 GeV 
and nib = 4:. 3 GeV just as in the MS scheme. 

That is, our Model consists of a set of equations (13) with (10) - (12) and 



= L2(a = 2) = L + 51nVZ^T2^, = ^2(2) = I + BlnVP + 27t^ . (13) 

Here, L and / contain common values, like in ( fTOl) . for each of the flavor region. 
Advantage of Model (fT0 l) -(|TT i) .(|T2l). ( |T3i) is that it involves only one new parameter, a = 2 
with Aj^ and nj taking their usual values. 

4.2 Accuracy of the Model vs data errors 

In paper \l\ errors of Model expressions (fTTl) . f|T2|) . f|T3l) in each Uf range were estimated 

by numerical comparison with the Magradze tables in the interval of 3-loop A^^"^^ ~ 
350 - 400 MeV . 

As it follows from that analysis, errors of Model for the first three APT functions are 
small, being of an order of 1-2 per cent for the first functions, of 3-5% for the second and 
of 6-10% for the third ones in the region above 1.5 GeV, i.e., in the Uf = 4,5 ranges. 
However, its accuracy in the Uf = 3 region (above IGeV ) is at the level of 5-10 per cent. 
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Meanwhile, relative contributions of typical LO, NLO and NNLO terms in APT non- 
power expansion for observables are usually something like 60-80%, 30-10%, and 10-1%, 
respectively (see Table 2 in Ref.[S]). Due to this, the Model accuracy for many cases is 
defined by that of the first model functions A^"'^ , W^"'^ , provided that QCD contribution 
to an observable starts from one-loop contribution ~ . At the same time, for quarko- 
nium decays the leading contribution ~ . There, the Model error is defined by accuracy 
of the third Minkowskian function Sl^"'^ . 

In the Table 2 of paper [T] we compared Model errors with some data errors in the 
low energy region. 

With due regard for data errors, we can now set some total margin of accuracy that 
our Model would satisfy. This margin may be chosen, e.g., as 1/3 of the data error bar, 
which is no less than 10%. Then, the accuracy limit, imposed on the first APT functions 
could be about 3%, for the second function ca 10%, and for the third ones, about 20%. 
Due to this, it is not reasonable to use the Model below 0.5 GeV , whereas above this 
limit it is fully advisable. Now, we proceed to its practical application to T decay. 



5 T(1S) Decay Revised 

5.1 T Widths 

Our first motive for this choice is that the parameter eA/(L) = vr^/L^ responsible for 
deviation of APT Minkowskian functions from powers of canonical as , is not very small in 
the region (5-10 GeV) related to this decay: eAf(L) ~ 0.16-0.27, £4/(^2) ~ 0.11-0.18. 

Another processes involving T , are T radiative decays and T production. But they 
have low data precision (about 10% — 40%). Besides, the non-radiative decay of IS-state 
provides the best data precision (1.5% for the ratio of hadronic and leptonic widths [3]). 

The second arguument is the disagreement with the world average. 

For an extensive review on qq decay widths see |T7]. NLO ratio of hadronic and 
leptonic decay widths of S state of the T was given in [18]. However, their expression is 
not renorm-invariant. In our analysis, to return it to the RG-invariant form, we put scale 
parameter ^ = My . Then 

r (T -> hadrons) 10(7r^ - 9)a^(gT) (, , a.(gY) ^ ^\ ,.2 (..^ 

1 (i^e+e ) 9 7ra^(Mx) V ^ / 

Then, the issue of scale should be readdressed to choice of st • 

5.2 Reevaluation of Aqcd from T Decay 

First attempt to evaluate Aqcd extracted from T decay by proper taking into account 
analytic continuation effects was made in [TT]. Analogous analysis is performed here. 
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employing better accuracy Model expression for 2tj(s) and using more fresh data from 
CLEO III detector [TH] as well as APT expansion instead of formula ([SD, that, within 
our Model, is reduced to 

5360 (4'^(L2) + 2.30 2li'^(L2)) . (15) 

We extracted by this formula and as values, from fresh data = 37.3±0.75|19j. 
In Table , results, obtained within APT are compared to results of standard PT. In line 

Table 1: Results of various as extraction from Upsilon decays 



Part I. Non-APT treatment 


Source 


a{Mr) 


aiMz) 


^%7^ 


PDG, T,1S 


0.170(4) 


0.112(2) 




PDG, global Fit 


0.182(5) 


0.1185(20) 


217+^g 


Part II. APT treatment 


Exact APT, IS 


0.1805(12)e.p 


0.1179(5)e.p 


210(5) 


[Mod], IS 


0.185(5)m 


0.120(2)m 


235(25) 


[5] Crude APT 


0.183 


0.119 


222 



1 of Part I, "PDG, IS"', standard PT results on T(IS') decay are given. We present them 
not exactly as they were published in [3] but recalculated along with modern experimental 
data. Line 2, marked "PDG, Fit" gives the published world average value described by 
the curve on Fig. 9.2 in [3J, within the error bars of all the processes. Column "a(Mx)" 
means "a^, calculated at the mass of T, according to eq.(9.5) of |3]". 

Line 1 of Part II, "Exact APT", presents results of TIS decay calculated by exact 
numeric tables for Sis and by [20] for 2I4 APT function. Line 2, [Mod], presents values 
obtained from T(IS') decay data by means of the Model eqs.(12),(13). Here, model errors 
combine Model errors of both the terms in the r.h.s. of eq.(15). Line 3, "Crude APT" 
gives an earlier result [5] with crude APT estimate used to correct the Bethke-2000 value 
as{Mr) = 0.170 extracted there from all the T decays data. 
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